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### / 

Abstract. This paper investigates the effect of collisional enhancement of the hydrogen lines on the derivation 
of the helium abundances in low metallicity H II regions. For this, we have constructed a grid of photoionization 
models relevant for the analysis of giant H II regions in blue compact galaxies. We show that the effect of 
collisional excitation on the Hq/H/3 ratio can be quite important (up to 8% or more). The impact of this effect 
on the determination of the helium mass fraction has been tracked on four low-metallicity blue compact galaxies 
for which Keck spectra are available and which are among the best objects for the quest of the pregalactic helium 
abundance. We find that taking into account the effects of collisional excitation of hydrogen results in an upward 
correction of the helium mass fraction Y by up to 5%. However, combining with other systematic effects usually 
not considered in the determination of the helium abundance in low-metallicity galaxies, the resulting uncertainty 
should be much less. 
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1. Introduction 

Tests of the big bang models and the determination of cos- 
mological parameters require an accurate determination of 
the relative abundances of the light isotopes produced dur- 
ing the big bang (Steigman et al. 1977; Walker et al. 1991; 
Sarkar 1996). The abundance of primordial 4 He is mostly 
derived from the helium abundances in metal-poor extra- 
galactic H II regions (Kunth & Sargent 1983; Pagel et al. 
1992; Izotov & Thuan 1998) although other nebulae such 
as H II regions in the Magellanic Clouds (Peimbert et al. 
2000) or planetary nebulae in the Galactic halo (Peimbert 
1983, 1989; Clegg 1989) have also occasionally been used. 
This is a challenging task, since the required accuracy is 
on the one percent level. 

In order to achieve such a goal, excellent quality ob- 
servations are obviously required together with an ex- 
tremely reliable reduction procedure. With the advent of 
very large telescopes and the progress in detectors, recent 
years have produced a large amount of data suitable for 
such an enterprise (Izotov et al. 1997; Izotov & Thuan 
1998; Izotov et al. 1999). As regards the interpretation of 
the observed line intensities in terms of abundance ratios, 
many problems occur and still remain to be overcome. 
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For example, the atomic data relevant for the helium line 
emissivities have only recently reached a state compati- 
ble with the demand (Benjamin et al. 1999). Among the 
problems to be solved for the determination of the he- 
lium abundance in one nebula, and towards which various 
groups of researchers have had different approaches, are: 
the stellar absorption lines underlying the nebular emis- 
sion lines, the reddening correction, the distribution of 
the electron temperature inside the nebulae, the density 
structure and its effect on the line emissivities, the optical 
depth in the lines, the ionization correction to derive the 
He/H ratio from the measured He + /H + ratio. The next 
step is to derive the pregalactic helium abundance from 
the measurement of the helium abundance in a sample of 
galaxies, by extrapolating the results towards zero metal- 
licity. These problems have been extensively discussed in a 
long list of papers including Peimbert & Torrcs-Pcimbcrt 
(1974, 1976), Stasinska (1980), Rayo et al. (1982), Kunth 
& Sargent (1983), Davidson & Kinman (1985), Dinerstein 
& Shields (1986), Pagel et al. (1986, 1992), Campbell 
(1992), Skillman & Kennicutt (1993), Balbes et al. (1993), 
Mathews et al. (1993), Olive & Steigman (1995), Sasselov 
& Goldwirth (1995), Olive et al. (1997), Steigman et al. 
(1997), Izotov & Thuan (1998), Izotov et al. (1997, 1999), 
Armour et al. (1999), Ballantyne et al. (2000), Viegas et 
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al. (2000), Peimbert et al. (2000, 2001), Olive & Skillman 
(2001), Sauer & Jedamzik (2001). 

The present day situation is that different groups reach 
different estimates of the pregalactic helium abundance, 
mutually exclusive within the quoted error bars. For ex- 
ample, Olive et al. (1997) give Y p = 0.234±0.002, while 
Izotov & Thuan (1998) give Y p = 0.244±0.002. These two 
estimates lead to very different cosmological implications, 
as discussed by Izotov et al. (1999). 

The effect on the pregalactic helium abundance deter- 
mination arising from deviation from case B theory for the 
hydrogen lines has been relatively little discussed so far. 
The effect of finite optical depth and dust on H/3 emis- 
sivity have been discussed by Cota & Ferland (1988) and 
Hummer & Storey (1992). The effect of collisional excita- 
tion of the hydrogen lines by thermal electrons has been 
pointed out by Davidson & Kinman (1985) and discussed 
for the case of I Zw 18 by Sasselov & Goldwirth (1995) and 
Stasinska & Schaerer (1999). Stasinska (2001) has shown 
that this effect can lead to an underestimation of the He 
I 5876/H/3 ratio by as much as 10% in hot planetary neb- 
ulae, independent of the actual reddening of the object. 

The purpose of the present paper is to systematically 
explore the effects of collisional excitation of the hydrogen 
lines on the helium abundance determinations in metal- 
poor giant H n regions. For this, we use a grid of photoion- 
ization models of giant H II regions and consider them in 
the light of the spectral properties of a sample of blue 
compact H n galaxies. This allows us to obtain a realistic 
estimate of the effect of collisional excitation by thermal 
electrons both on Ha and on H/3. We then concentrate 
on a selection of blue compact galaxies with particularly 
good quality spectroscopic data, and track the effect of 
collisional excitation - at the level predicted by the mod- 
els - on the derived helium abundance for these objects. 

2. The photoionization models 

The photoionization models of giant H II regions are con- 
structed using the procedure described in Stasinska et al. 
(2001, hereinafter SSL2001) to which we refer for details. 
Briefly, we consider coeval stellar clusters of given ini- 
tial mass M+ that ionize the surrounding gas, supposed 
to be of same chemical composition as the stars. We as- 
sume instantaneous bursts with a Salpeter initial mass 
function and an upper stellar limit M up = 120 M©. The 
radiation field is provided by the evolutionary synthesis 
models of Schaerer & Vacca (1998) based on the non- 
rotating Geneva stellar evolution models, with the high 
mass-loss tracks of Meynet et al. (1994). The spectral 
energy distributions for massive main-sequence stars are 
those given by the CoStar models (Schaerer & de Koter 
1997). The pure He models of Schmutz ct al. (1992) are 
used for Wolf-Rayet stars. The ionization and temperature 
structure of the nebula is computed with the photoioniza- 
tion code PHOTO, as described in Stasinska & Leitherer 
(1996, hereinafter SL96). We assume that the nebulae are 
spherical and radiation-bounded, with uniform density n 



= 10 cm" 3 and filling factor e — 1 and that all the stars 
are located at the center. As in our previous papers (SL96 
and SSL2001) and in all the papers based on static pho- 
toionization models, our models do not consider the pos- 
sible pollution of the nebulae by stellar ejecta during the 
course of the evolution of the star clusters. While this 
phenomenon may have important consequences, it is not 
really relevant to the problem under study here. 

The grid used in the present study is composed of 
models with metallicities Z of 0.2, 0.05 and 0.02 times 
the solar metallicity (using the prescription of McGaugh 
(1991) for the relative abundances of the elements). As 
is known (e.g. Shields 1986) the emission line spectrum 
of photoionized nebulae is essentially dependent on the 
spectral energy distribution of the ionizing radiation field, 
on the metallicity and on the ionization parameter U. For 
a constant density sphere U is equal to A(Q(H°)ne 2 ) 1 / 3 , 
where Q(H°) is the total number of H Lyman continuum 
photons emitted by the ionizing source (proportional to 
the mass of the star cluster) and A is approximately given 
by A w 2.8 x 10- 20 (10 4 /T C ) 2/3 , where T c is the electron 
temperature. At a given chemical composition, any com- 
bination of Q(H ), n and e giving the same value of U will 
result in the same emission line spectrum (at least when 
the densities are sufficiently low that cooling is not affected 
by collisional deexcitation of forbidden lines, which is the 
case of the objects under study in this paper). Three dif- 
ferent initial masses for the star clusters are considered: 
M*=1Q 3 , 10 6 , and 10 9 M Q (assuming a lower IMF mass 
cut-off of Mi ow = 0.8M Q ). With such a grid, we reasonably 
cover the range of metallicities and ionization parameters 
relevant to the objects we are interested in. 

For each model, we derive the electron tempera- 
ture T[ IHJ4363/5007 from the computed [O iii]A4363/[0 
ni]A5007 ratioR, using the same atomic data as in the 
photoionization code. Figure 1 presents the sequences 
of models in the Tjo m]4363/5007 versus [O iii]A5007/[O 
n]A3727 plane (left) and in the Tjo ni]4363/5007 versus [O 
iii]A5007/H/3 plane (right). The top panels correspond 
to a total initial stellar mass of M*=1O 3 M0, the middle 
ones to Af*=10 6 M Q and the bottom ones to M*=1O 9 M . 
Symbols mark the epoch by time steps of 10 6 yr, starting 
from an age of 10 4 yr. Series with Z= 0.02 Zq are repre- 
sented by circles, series with Z= 0.05 Zq by squares, and 
series with Z= 0.2 Zq by triangles. The time evolution is 
from the right to the left in these plots, since the ionization 
parameter decreases due to gradual disappearance of the 
most massive stars. This provokes a softening of the ion- 
izing radiation field, which induces a decrease in the elec- 
tron temperature. As expected, the highest temperatures 
are found in the models with lowest metallicities (cooling 
by collisional excitation of forbidden lines becomes less 
important) and highest ionization parameters (cooling by 
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of the line at 5007A, not to the sum of the lines [O iii]A4959 
and A5007. 
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Fig. 1. Sequences of photoionization models for instantaneous bursts of various metallicities. The metallicity is indi- 
cated by the following symbols: circle for Z /Zq = 0.02, square for Z / 'Zq = 0.05, triangle for Z j Zq = 0.2. The symbols 
mark time steps of 1 Myr. The models shown in the top panels have M*= 10 3 M , those in the middle panels have 
M*= 10 6 Mq, those in the bottom panels have M*= 10 9 Mq (see text). The ordinate represents the temperature 
indicated by the [O iii]A4363/[0 iii]A5007 line ratio. In the left panels it is shown as a function of [O iii]A5007/[O 
n]A3727, in the right panels as a function of [O m]A5007/H/3. 



collisional excitation of H Lyman a becomes less impor- 
tant). 

In the next section, we discuss the relevance of this 
grid of models for the problem we are investigating in the 
present paper. 

3. The sample of blue compact galaxies 
confronted to the model grid 

Since ab initio models of nebulae photoionized by evolv- 
ing starbursts do not necessarily well represent the ob- 



served properties of giant H n regions (see discussion in 
SSL2001), we have to use them in the light of observed 
properties of H n galaxies. For this, we consider a sample 
of blue compact galaxies observed since 1994 by Izotov 
and coworkers. In addition to the 69 objects considered in 
the study by SSL2001, it contains 64 new objects observed 
with different telescopes, including Keck, VLT, MMT, 4m 
and 2.1m KPNO (Izotov et al. 2001a,b; Guseva et al. 2001; 
Papaderos et al. 2001). The new observations were re- 
duced in the same way as those considered by SSL2001. 
In all cases, the line intensities have been corrected for 
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Fig. 2. Our observational sample in various diagrams relating emission line ratios and H/3 equivalent width. On the 
upper right of each panel is given the total number of the objects appearing in the diagram. The behaviour of the 
present sample in these diagrams is very similar to the behaviour of the smaller sample shown by Stasihska et al. 
(2001) in their Fig. 2. 



reddening by using a procedure described in Izotov et al. 
(1994) in which the reddening constant and stellar absorp- 
tion in the Balmcr lines are determined simultaneously. 
Figure 2 shows the same observational diagrams as Fig. 
2 of SSL2001 for the new sample of 133 objects. Exactly 
the same trends are seen (with a slightly higher dispersion 
due to the fact that some of the new data are of lower 
quality than those of the sample considered by SSL2001). 



Therefore, we are dealing with the same population of ob- 
jects. Hence the conclusions of SSL2001 still apply: the 
sample is mostly composed of objects whose ionization is 
produced by a starburst of age younger than about 5 Myr. 
In the following, we will then restrict our considerations 
to models corresponding to such ages. 

Figure 3 shows the observed sample in the 
T[o ni]4363/5007 versus [O iii]A5007/[O n]A3727 plane (left) 
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Fig. 3. Our observational sample in the same planes as the models shown in Fig. 1. 
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and in the Tjo m]4363/5007 versus [O iii]A5007/H/3 plane 
(right). By comparing with Fig. 1, we see that while the 
models with the highest Tr 111)4363/5007 wen reproduce the 
[O m]A5007/H/3 ratios, they have [O in]A5007/[O n]A3727 
ratios much larger than observed in objects of similar 
values of 7] 111)4363/5007 • For the range of observed val- 
ues of [O iii]A5007/[O n]A3727, it is impossible to obtain 
models with 7] 111)4363/5007 around 20 000K. One possi- 
ble way to increase the [O 11] line intensity is to consider 
composite models, in which part of the emission comes 
from higher density filaments. However, it was impossi- 
ble to find a prescription for the composite models to re- 
produce at the same time [O m]A4363/[0 m]A5007, [O 
lll]A5007/[O n]A3727 and the electron density sensitive 
ratio [S n]A6731/6717. The conclusion is that simple pho- 
toionization models are unable to reproduce the high [O 
iii]A4363/[0 m]A5007 ratios that are observed in quite a 
number of giant H II regions. A similar conclusion has 
been reached from tailored model fitting of various gi- 
ant H II regions (Campbell 1988; Stasihska & Schaerer 
1999; Luridiana et al. 1999; Luridiana & Peimbert 2001). 
Reducing the abundances of coolants such as carbon or 
silicon (on which there are virtually no observational con- 
strains) is insufficient to solve the problem. It seems that 
an additional heating mechanism is at work, and a popu- 
lar candidate is shock heating. Whatever may be the ori- 
gin of the additional heating, this means that simple pho- 
toionization models returning a [O iii]A5007/[O n]A3727 
in the observed range are likely to have too low an electron 
temperature, and that those returning a [O iii]A4363/[0 
iii]A5007 comparable with observations have too high a 
degree of ionization. This means that the amount of col- 
lisional excitation they predict for the hydrogen lines is 
likely a lower limit, since the contribution of collisional 
excitation increases with electron temperature and with 
the proportion of residual neutral hydrogen. 



4. Model predictions for the collisional excitation 
of hydrogen lines 

In Figure 4, we display the same sequences of models as in 
Figure 1, but up to ages of 5 Myr, in the T [Q m] 4 363/5007 
versus [O lii]A5007/H/3 plane. Again, the top panels cor- 
respond to a total initial stellar mass of M*=1O 3 M , 
the middle ones to M^—10 6 M Q and the bottom ones to 
M*=10 9 Mq. The collisional contribution to the Balmcr 
line emission is indicated in the following way. In the left 
panels, next to the models are indicated the values of : 



1000 x 



/ (Ha/H/3) (model) 
^(Ha/H/3)(casc B) 



where the case B value of Ha/H/3 is computed for the tem- 
perature Tjo iii]4363/5007 deduced from the [O 111] A4363/[0 
iii]A5007 ratio predicted by the model. 

It is seen that in our models the Ha/H/3 ratios may 
be larger than the case B value (at the temperature 
T[o HIJ4363/5007) by up to 8%. This means that the red- 
dening computed assuming the case B value of the Ha/H/3 
ratios may be significantly overestimated. 

In the right panels, next to the models, we indicate 
(in units of 1 per 1000) the collisional contribution to the 
emission in the H/3 line: 



1000 x 



/H/3(modcl) 
l,H/3(case B) 



1 



where H/3(case B) stands for the integrated emission of the 
model in the H/3 line at the local electron temperature. 

Clearly, the collisional contribution to H/3 in our mod- 
els may reach 2%. 

The computations in PHOTO have been made with 
the collisional excitation rate coefficients from Aggarwal 
(1983) for the n — 3 level of hydrogen, and from Drake & 
Ulrich (1980) for higher levels. Values given by other au- 
thors can be significantly different. For example Aggarwal 
et al. (1991) give values that are larger than the ones from 
Drake & Ulrich by up to a factor 2. Unfortunately, there is 
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Fig. 4. The same models as in Fig. 1 shown in the Xjo m]4363/5007 versus [O iii]A5007/H/3 plane, but only for starburst 
ages up to 5 Myr. As in Fig. 1, the models shown in the top panels have M*= 10 3 M Q , those in the middle panels 
have M*= 10 6 M , those in the bottom panels have M±= 10 9 Mq. The numbers indicate, in units of 1 per 1000, the 
values of deviation from case B theory for Ha/H/3 (left panels) and B./3 (right panels). See text for a more detailed 
explanation. 



so far no agreement on what are the most reliable values 
of the collision strengths for hydrogen (Callaway 1994). 

Note that there is no simple relation between the to- 
tal Ha and H/3 emission in our models, although at each 
point in the nebula the relation between the Ha and H/3 
emissivities is a function of the local temperature only. 
The reason is that the electron temperature in our models 
is not uniform, and the relation between local Ha and H/3 
emissivities do not propagate in a simple way into inte- 
grated line ratios. It is for a similar reason that a temper- 
ature correction factor tcf has been introduced by Sauer 



& Jedamzik (2001) to derive the helium abundance from 
the ratios of helium to hydrogen line intensities. 

In real nebulae, where the temperature structure is 
much more erratic than in simple photoionization models, 
the situation is probably even more complicated. 

One might be tempted to use the ratio of H Lyman 
a to H/3 to better constrain the amount of collisional ex- 
citation of hydrogen lines in real nebulae. However, this 
proves difficult, because the path of H Lyman a photons 
is strongly increased by scattering in the ionized as well as 
in the neutral part of the nebula, during which absorption 
by dust and molecular hydrogen occur. Therefore, the ob- 
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served intensity of H Lyman a strongly depends on the 
gas and dust distribution and on the velocity field. 

5. The effect of H collisional excitation on the 
derived helium abundances for a few selected 
blue compact galaxies 

It is instructive to investigate what are the effects of col- 
lisional excitation of the hydrogen lines on the derived 
helium abundance. This is done in the present section on 
a few examples. 

We have selected 4 low-metallicity blue compact galax- 
ies SBS 0335-052 (Z©/40, Izotov et al. 2001b), SBS 
0940+544 (Z©/27, Guseva et al. 2001), Tol 1214-277 and 
Tol 65 (Zq/23, Izotov et al. 2001a) from the sample pre- 
sented in Section 3, on the basis of their high-quality spec- 
tra obtained with the Keck telescope and high electron 
temperatures in the range of 17000 - 20000 K derived from 
the [O in] line intensity ratios. Because of the high electron 
temperature the effect of collisional excitation on the he- 
lium abundance determination is expected to be relatively 
large in these galaxies. To derive the 4 He and heavy ele- 
ment abundances we follow the prescriptions by Izotov et 
al. (1994), but assuming that there is a certain percentage 
of collisional excitation in the Ha and in the H/3 lines. We 
assume that the temperature characteristic of the emis- 
sion of all He I lines is given by Tiq iii]4363/5007j an d that 
the ionization correction factor for helium is equal to 1. 
This is not exactly true for all our models, as seen in Fig. 
5, which shows (Tjo ni]4363/5007/^ 1 (H + ) (left panels) and 
(.t(Hc+) + .t(Hc ++ ))/:e(H + ) (right panels) for the same 
models as Fig. 4. T(H + ) is defined as 



T e (H+) = 



/ T e n(U+)n e dV 
Jn(R+)n e dV 



and x(Hc + ) is defined as 



/ n(Rc + )n e dV 



Jn(Re)n e dV ' 



(1) 



(2) 



However, since the four considered galaxies have 
T[o iii]4363/5007 around 20000 K and [O m]A5007/H/3 
larger than 3, our model grid indicates that our approxi- 
mations should be reasonable. 

In order to explore the effects of collisional excitation 
of the H lines, we consider two different cases. In the first 
case the effect of collisional excitation is accounted for by 
decreasing by some amount the intensity of the Ha emis- 
sion line and leaving unchanged the intensity of the H/3 
emission line. In the second case we decrease the intensities 
of Ha and H/3 emission lines assuming that the effect of the 
collisional excitation on the H/3 line intensity is three times 
smaller than that on the Ha line intensity. For each as- 
sumed value of the contribution of the collisional enhance- 
ment to Ha, A/(Ha)//(Ha), we compute the extinction, 
correct the emission line ratios for reddening and under- 
lying stellar absorption, then recompute TJq iii]4363/5007 
and derive the 4 He and the heavy element abundances. In 



the analysis, we use the five strongest He I emission lines 
to derive the helium mass fraction Y, after correction for 
collisional and fluorescent excitation of the He I emission 
lines. In short, using the intensities of the five He I emission 
lines, we evaluate the electron number density n e in the 
He + zone and the optical depth t(A3889) in the He I A3889 
line in a self-consistent way, so that the He I A3889/A5876, 
A4471/A5876, A6678/A5876 and A7065/A5876 line ratios 
have their recombination values, after correction for col- 
lisional and fluorescent enhancement of the He I emission 
lines. These effects are very strong in hot low-metallicity 
H II regions. In particular, the correction of the He I A5876 
emission line for collisional and fluorescent enhancement 
in SBS 0335-052 is - 12%. 

Of course, for each object, there is a maximum con- 
tribution of collisional enhancement A/(Ha)//(Ha) im- 
posed by the observed Ha/H/3. In reality, the range of 
possible values of A7(Ha)//(Ha) is smaller because in- 
terstellar extinction must be present at some level. We 
have checked that, after the completion of the whole pro- 
cedure, the higher order Balmer lines remain compatible 
with case B within observational errors and uncertainties 
in the extinction law. 

In Fig. 6 we show the variations of the correction of 
the helium mass fraction AY/Y derived from the He I 
A5876 emission line for the four galaxies we investigated. 
It is seen that the correction for the collisional enhance- 
ment of the hydrogen lines results in an increase of the 
He abundance in all galaxies. Qualitatively such a depen- 
dence can be explained by two effects: 1) the correction 
of the H/3 line in the second case (dotted line) increases 
the He i/H/3 intensity ratios; 2) the correction of Ha and 
H/3 lines for collisional excitation decreases the extinction 
and hence increases the intensity of He I A5876 emission 
line. Certainly, some part of the deviation of the observed 
Ha/H/3 intensity ratio from the case B value is due to 
interstellar extinction. Therefore, although the effect of 
collisional excitation of the hydrogen lines on the helium 
abundance is significant, it is unlikely that the correction 
AY/Y exceeds 5%. At larger values of AY /Y the inter- 
stellar extinction derived from the Ha/H/3 becomes very 
small, less than the foreground extinction produced in our 
Galaxy. Our estimate of AY/Y is not very certain and 
should be considered as a qualitative indication of the rel- 
atively large positive correction of the helium abundance 
if collisional excitation of the hydrogen emission lines is 
taken into account. An accurate determination of the he- 
lium abundance requires detailed photoionization model- 
ing of the H II regions in individual galaxies, with all rel- 
evant effects taken into account, and with all the relevant 
line intensities satisfactorily reproduced. Such an enter- 
prise is still ahead of us. 

6. Conclusion and prospects 

Our consideration of the collisional excitation of hydrogen 
emission lines has shown that this effect can be particu- 
larly large, as high as 8% for the Ha/H/3 ratio in the low- 
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Fig. 5. The same models as in Fig. 4. The symbols have the same meaning as in Fig. 1. As in Figs. 1 and 4, models 
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helium and show (x(Hc + ) + .t(H + ))/x(H+) versus [O lll]A5007/H/3. 



metallicity extragalactic H n regions. If taken into account 
it results in an upward correction of the helium abundance 
Y by up to 5%, making the collisional excitation of hy- 
drogen one of the most important sources of systematics 
in the primordial helium abundance determination. Some 
other sources of systematic uncertainties should be consid- 
ered as well. One of the most important sources may be 
the "temperature fluctuations" which , if present, result in 
the downward correction of Y by ^ 5% (e.g., Peimbert et 
al. 2000). Sauer & Jedamzik (2001) considered the effect of 
the ionization structure and large-scale temperature vari- 
ations in the H n region on the He abundance determi- 
nation. Izotov et al. (2001a) compared the observed prop- 



erties of the high-excitation H II regions in the galaxies 
Tol 1214-277 and Tol 65 with Sauer & Jedamzik models 
and found that a small downward Y correction (by less 
than 1%) should be applied due to the two mentioned ef- 
fects. Finally, underlying stellar He I line absorption of the 
ionizing clusters should be taken into account. However, 
this effect is not large in the four considered blue compact 
galaxies because of the large equivalent widths of emis- 
sion lines and results in an upward Y correction by S_. 
1%. Combining all these systematic effects it appears that 
the helium abundance derived by Izotov & Thuan (1998) 
and Izotov et al. (1999) in low-metallicity blue compact 
galaxies could be underestimated by about 2 - 3% only. 
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However, additional work has to be done to clarify the im- 
portance of the systematic effects in the 4 He abundance 
determination. Especially, the mystery of the "tempera- 
ture fluctuations" has to be solved. For this, a detailed 
modeling of the H n regions in each of the best observed 
low-metallicity blue compact galaxies should be done. 
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